Genomic and Phenotypic Analysis of Salmonella enterica Bacteriophages Identifies Two Novel Phage Species.

Bacteriophages (phages) are potential alternatives to chemical antimicrobials against pathogens of public health significance. Understanding the diversity and host specificity of phages is important for developing effective phage biocontrol approaches. Here, we assessed the host range, morphology, and genetic diversity of eight Salmonella enterica phages isolated from a wastewater treatment plant. The host range analysis revealed that six out of eight phages lysed more than 81% of the 43 Salmonella enterica isolates tested. The genomic sequences of all phages were determined. Whole-genome sequencing (WGS) data revealed that phage genome sizes ranged from 41 to 114 kb, with GC contents between 39.9 and 50.0%. Two of the phages SB13 and SB28 represent new species, Epseptimavirus SB13 and genera Macdonaldcampvirus, respectively, as designated by the International Committee for the Taxonomy of Viruses (ICTV) using genome-based taxonomic classification. One phage (SB18) belonged to the Myoviridae morphotype while the remaining phages belonged to the Siphoviridae morphotype. The gene content analyses showed that none of the phages possessed virulence, toxin, antibiotic resistance, type I-VI toxin-antitoxin modules, or lysogeny genes. Three (SB3, SB15, and SB18) out of the eight phages possessed tailspike proteins. Whole-genome-based phylogeny of the eight phages with their 113 homologs revealed three clusters A, B, and C and seven subclusters (A1, A2, A3, B1, B2, C1, and C2). While cluster C1 phages were predominantly isolated from animal sources, cluster B contained phages from both wastewater and animal sources. The broad host range of these phages highlights their potential use for controlling the presence of S. enterica in foods.

Survival and Expression of rpoS and grxB of Cronobacter sakazakii in Powdered Infant Formula Under Simulated Gastric Conditions of Newborns.

Cronobacter sakazakii can cause severe illnesses in infants, predominantly in preterm newborns, with consumption of contaminated powdered infant formula (PIF) being the major vehicle of infection. Using a dynamic human gastrointestinal simulator called the SHIME, this study examined the effects of gastric acidity and gastric digestion time of newborns on the survival and expression of stress genes of C. sakazakii. Individual strains, inoculated at 7 log CFU/mL into reconstituted PIF, were exposed to gastric pH values of 4.00, 5.00 and 6.00 for 4 h with gradual acidification. The survival results showed that C. sakazakii grew in the stomach portion of the SHIME during a 4-h exposure to pH 4.00, 5.00 and 6.00 by 0.96-1.05, 1.02-1.28 and 1.11-1.73 log CFU/mL, respectively. The expression of two stress genes, rpoS and grxB, throughout gastric digestion was evaluated using reverse transcription qPCR. The upregulation of rpoS and grxB during the 4-h exposure to simulated gastric fluid at pH 4.00 showed that C. sakazakii strains may be experiencing the most stress in the pH 4.00 treatment. The gene expression results also suggest that C. sakazakii strains appeared to develop an acid adaptation response during the 4-h exposure that may facilitate their survival. Altogether, this study highlights that a combination of low gastric acidity, long digestion time in the presence of reconstituted PIF, created a favorable environment for the adaptation and survival of C. sakazakii in the simulation of a newborn's stomach. This study gives directions for future research to further advance our understanding of the behavior of C. sakazakii in the GI tract of newborns.

Real-time evaluation of signal accuracy in wastewater surveillance of pathogens with high rates of mutation.

Wastewater surveillance of coronavirus disease 2019 (COVID-19) commonly applies reverse transcription-quantitative polymerase chain reaction (RT-qPCR) to quantify severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA concentrations in wastewater over time. In most applications worldwide, maximal sensitivity and specificity of RT-qPCR has been achieved, in part, by monitoring two or more genomic loci of SARS-CoV-2. In Ontario, Canada, the provincial Wastewater Surveillance Initiative reports the average copies of the CDC N1 and N2 loci normalized to the fecal biomarker pepper mild mottle virus. In November 2021, the emergence of the Omicron variant of concern, harboring a C28311T mutation within the CDC N1 probe region, challenged the accuracy of the consensus between the RT-qPCR measurements of the N1 and N2 loci of SARS-CoV-2. In this study, we developed and applied a novel real-time dual loci quality assurance and control framework based on the relative difference between the loci measurements to the City of Ottawa dataset to identify a loss of sensitivity of the N1 assay in the period from July 10, 2022 to January 31, 2023. Further analysis via sequencing and allele-specific RT-qPCR revealed a high proportion of mutations C28312T and A28330G during the study period, both in the City of Ottawa and across the province. It is hypothesized that nucleotide mutations in the probe region, especially A28330G, led to inefficient annealing, resulting in reduction in sensitivity and accuracy of the N1 assay. This study highlights the importance of implementing quality assurance and control criteria to continually evaluate, in near real-time, the accuracy of the signal produced in wastewater surveillance applications that rely on detection of pathogens whose genomes undergo high rates of mutation.

Whole genome sequence of Vibrio cholerae NB-183 isolated from freshwater in Ontario, Canada harbors a unique gene repertoire.

OBJECTIVE: Vibrio cholerae is an enteric pathogen that poses a significant threat to global health. It causes severe dehydrating diarrheal disease cholera in humans. V. cholerae could be acquired either from consuming contaminated seafood or direct contact with polluted waters. As part of a larger program that assesses the microbial community profile in aquatic systems, V. cholerae strain NB-183 was isolated and characterized using a combination of culture- and whole-genome sequencing-based approaches. DATA DESCRIPTION: Here we report the assembled and annotated whole-genome sequence of a V. cholerae strain NB-183 isolated from a recreational freshwater lake in Ontario, Canada. The genome was sequenced using short-read Illumina systems. The whole-genome sequencing yielded 4,112,549 bp genome size with 99 contigs with an average genome coverage of 96× and 47.42% G + C content. The whole genome-based comparison, phylogenomic and gene repertoire indicates that this strain harbors multiple virulence genes and biosynthetic gene clusters. This genome sequence and its associated datasets provided in this study will be an indispensable resource to enhance the understanding of the functional, ecological, and evolutionary dynamics of V. cholerae.

Draft genome sequences of two Salmonella Uzaramo isolates from poultry in South Africa.

Salmonella enterica is a zoonotic pathogen and a leading cause of foodborne gastroenteritis in humans. Here, we report the draft genome sequences of two Salmonella Uzaramo isolates, which were isolated from poultry organs during routine post-mortem examination in South Africa. Currently, whole-genome sequences on Salmonella Uzaramo are scanty.

Draft genome sequences of two Proteus mirabilis isolates recovered from a municipal wastewater treatment plant in Ontario, Canada.

Proteus mirabilis is a Gram-negative bacterium that is frequently implicated in urinary tract infections in humans and companion animals and has also been associated with foodborne infections in several countries. Here, we report the draft genome sequences of two P. mirabilis isolates recovered from municipal wastewater.

Genetic diversity of Salmonella enterica isolated over 13 years from raw California almonds and from an almond orchard.

A comparative genomic analysis was conducted for 171 Salmonella isolates recovered from raw inshell almonds and raw almond kernels between 2001 and 2013 and for 30 Salmonella Enteritidis phage type (PT) 30 isolates recovered between 2001 and 2006 from a 2001 salmonellosis outbreak-associated almond orchard. Whole genome sequencing was used to measure the genetic distance among isolates by single nucleotide polymorphism (SNP) analyses and to predict the presence of plasmid DNA and of antimicrobial resistance (AMR) and virulence genes. Isolates were classified by serovars with Parsnp, a fast core-genome multi aligner, before being analyzed with the CFSAN SNP Pipeline (U.S. Food and Drug Administration Center for Food Safety and Applied Nutrition). Genetically similar (≤18 SNPs) Salmonella isolates were identified among several serovars isolated years apart. Almond isolates of Salmonella Montevideo (2001 to 2013) and Salmonella Newport (2003 to 2010) differed by ≤9 SNPs. Salmonella Enteritidis PT 30 isolated between 2001 and 2013 from survey, orchard, outbreak, and clinical samples differed by ≤18 SNPs. One to seven plasmids were found in 106 (62%) of the Salmonella isolates. Of the 27 plasmid families that were identified, IncFII and IncFIB plasmids were the most predominant. AMR genes were identified in 16 (9%) of the survey isolates and were plasmid encoded in 11 of 16 cases; 12 isolates (7%) had putative resistance to at least one antibiotic in three or more drug classes. A total of 303 virulence genes were detected among the assembled genomes; a plasmid that harbored a combination of pef, rck, and spv virulence genes was identified in 23% of the isolates. These data provide evidence of long-term survival (years) of Salmonella in agricultural environments.

Complete genome sequences of agricultural azole-resistant Penicillium rubens encoding CYP51A and ERG11 paralogues.

Azoles are major antifungals in agriculture and medicine. However, the surge of intrinsic azole resistance is critical for public health. Here, we present the complete long-read sequencing of three azole-resistant Penicillium rubens from food crops. The presence of CYP51A and ERG11 paralogues was confirmed, as in other azole-resistant P. rubens.

Evaluation of four human-associated fecal biomarkers in wastewater in Southern Ontario.

Human fecal biomarkers (HFBs) have a longstanding history in the field of microbial source tracking (MST) serving as indicators of human fecal contamination in drinking and recreational water. Further, HFBs have aided in recent efforts to monitor human pathogen transmission within communities. The dilution of wastewater from various sources throughout the sewershed cannot be controlled and human fecal biomarkers (HFBs) can be used to normalize target human pathogen concentrations so that fluctuations in fecal matter in wastewater can be accounted for. In the current study, we monitored the prevalence of four HFBs - including two viruses, Pepper mild mottle virus (PMMoV), cross-assembly phage (crAssphage), as well as two human-associated Bacteroides markers, HF183 and BacHuman - in wastewater samples from ten Southern Ontario wastewater treatment plants and evaluated their temporal and spatial variation in context of environmental factors that may impact the ability of HFB to normalize pathogen concentrations in wastewater. Environmental variables including precipitation, wastewater flow rate, temperature, and concentrated mass were also analyzed for their potential correlation with HFB variation in wastewater. The four HFBs were detected at high concentrations across all 10 sampling locations. The median concentrations across all sampling sites were: PMMoV 3.6 Log gene copies (GC)/mL; crAssphage 5.0 Log GC/mL; HF183 6.8 Log GC/mL and BacHuman 6.9 Log GC/mL. All HFBs were found to be similarly stratified across all 10 sites, and the bacterial markers were consistently found at higher concentration compared to the viral HFBs at all sites. The coefficient of variation (CV) for each HFB was used to characterize the variability of each biomarker at each sewershed. BacHuman and crAssphage were found to have lower CV than PMMoV and HF183, indicating that BacHuman and crAssphage may perform better in reflecting the variations in abundance of human feces in wastewater or MST applications.

Draft Genome Sequence of Exiguobacterium sp. Strain N5, Isolated from a Recreational Freshwater Kettle Lake in Ontario.

Exiguobacterium spp. are facultative anaerobic, Gram-positive, non-spore-forming bacilli, reported to tolerate extreme environments. Here, we report the draft genome sequence of Exiguobacterium sp. strain N5, isolated from a recreational freshwater lake.

Application of prophage sequence analysis to investigate a disease outbreak involving Salmonella Adjame, a rare serovar and implications for the population structure.

Introduction: Outbreak investigation of foodborne salmonellosis is hindered when the food source is contaminated by multiple strains of Salmonella, creating difficulties matching an incriminated organism recovered from patients with the specific strain in the suspect food. An outbreak of the rare Salmonella Adjame was caused by multiple strains of the organism as revealed by single-nucleotide polymorphism (SNP) variation. The use of highly discriminatory prophage analysis to characterize strains of Salmonella should enable a more precise strain characterization and aid the investigation of foodborne salmonellosis. Methods: We have carried out genomic analysis of S. Adjame strains recovered during the course of a recent outbreak and compared them with other strains of the organism (n = 38 strains), using SNPs to evaluate strain differences present in the core genome, and prophage sequence typing (PST) to evaluate the accessory genome. Phylogenetic analyses were performed using both total prophage content and conserved prophages. Results: The PST analysis of the S. Adjame isolates showed a high degree of strain heterogeneity. We observed small clusters made up of 2-6 isolates (n = 27) and singletons (n = 11) in stark contrast with the three clusters observed by SNP analysis. In total, we detected 24 prophages of which only four were highly prevalent, namely: Entero_p88 (36/38 strains), Salmon_SEN34 (35/38 strains), Burkho_phiE255 (33/38 strains) and Edward_GF (28/38 strains). Despite the marked strain diversity seen with prophage analysis, the distribution of the four most common prophages matched the clustering observed using core genome. Discussion: Mutations in the core and accessory genomes of S. Adjame have shed light on the evolutionary relationships among the Adjame strains and demonstrated a convergence of the variations observed in both fractions of the genome. We conclude that core and accessory genomes analyses should be adopted in foodborne bacteria outbreak investigations to provide a more accurate strain description and facilitate reliable matching of isolates from patients and incriminated food sources. The outcomes should translate to a better understanding of the microbial population structure and an 46 improved source attribution in foodborne illnesses.

Draft Genome Sequence of Bacillus anthracis N1, Isolated from a Recreational Freshwater Kettle Lake in Ontario, Canada.

Bacillus anthracis is widespread in soil and a causative agent of anthrax, primarily in herbivores. Here, we report the draft genome sequence of Bacillus anthracis strain N1, which was isolated from a recreational freshwater lake and found to carry multiple antibiotic resistance genes and biosynthetic gene clusters.

Inhibition of Cronobacter sakazakii in an infant simulator of the human intestinal microbial ecosystem using a potential synbiotic.

Powdered infant formula (PIF) can be contaminated with Cronobacter sakazakii, which can cause severe illnesses in infants. Synbiotics, a combination of probiotics and prebiotics, could act as an alternative control measure for C. sakazakii contamination in PIF and within the infant gut, but synbiotics have not been well studied for their ability to inhibit C. sakazakii. Using a Simulator of the Human Intestinal Microbial Ecosystem (SHIME®) inoculated with infant fecal matter, we demonstrated that a potential synbiotic, consisting of six lactic acid bacteria (LAB) strains and Vivinal GOS, can inhibit the growth of C. sakazakii in an infant possibly through either the production of antimicrobial metabolites like acetate, increasing species diversity within the SHIME compartments to compete for nutrients or a combination of mechanisms. Using a triple SHIME set-up, i.e., three identical SHIME compartments, the first SHIME (SHIME 1) was designated as the control SHIME in the absence of a treatment, whereas SHIME 2 and 3 were the treated SHIME over 2, 1-week treatment periods. The addition of the potential synbiotic (LAB + VGOS) resulted in a significant decrease in C. sakazakii levels within 1 week (p < 0.05), but in the absence of a treatment the significant decline took 2 weeks (p < 0.05), and the LAB treatment did not decrease C. sakazakii levels (p ≥ 0.05). The principal component analysis showed a distinction between metabolomic profiles for the control and LAB treatment, but similar profiles for the LAB + VGOS treatment. The addition of the potential synbiotic (LAB + VGOS) in the first treatment period slightly increased species diversity (p ≥ 0.05) compared to the control and LAB, which may have had an effect on the survival of C. sakazakii throughout the treatment period. Our results also revealed that the relative abundance of Bifidobacterium was negatively correlated with Cronobacter when no treatments were added (ρ = -0.96; p < 0.05). These findings suggest that C. sakazakii could be inhibited by the native gut microbiota, and inhibition can be accelerated by the potential synbiotic treatment.

Draft Genome Sequences of Two Clostridium botulinum Group II Strains Carrying Phage-Like Plasmids.

Clostridium botulinum is responsible for botulism, a potentially lethal foodborne intoxication. Here, we report the draft genome sequences of C. botulinum group II strains 202F (serotype F) and Hazen (serotype E). The genomes share many similarities, including multiple mobile genetic elements.

Metagenomics of Wastewater Influent from Wastewater Treatment Facilities across Ontario in the Era of Emerging SARS-CoV-2 Variants of Concern.

We report metagenomic sequencing analyses of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA in composite wastewater influent from 10 regions in Ontario, Canada, during the transition between Delta and Omicron variants of concern. The Delta and Omicron BA.1/BA.1.1 and BA.2-defining mutations occurring in various frequencies were reported in the consensus and subconsensus sequences of the composite samples.

Selection of a Potential Synbiotic against Cronobacter sakazakii.

ABSTRACT: Cronobacter sakazakii is an opportunistic foodborne pathogen that can be fatal to infants; it is commonly associated with powdered infant formula due to contamination during manufacturing processes or during preparation in hospitals or homes. This project aimed to select a potential synbiotic, a combination of probiotic strains with a prebiotic product, to inhibit the growth of C. sakazakii in an in vitro dynamic infant gut model (Simulator of the Human Intestinal Microbial Ecosystem). A total of 16 lactic acid bacteria (LAB) were tested for their inhibitory properties against four different C. sakazakii strains by a zone of inhibition test. Lactobacillus and Pediococcus species were able to inhibit the growth (>15-mm inhibition zones) of all C. sakazakii strains tested, and only one strain from the two genera exhibited atypical resistance to tetracycline. All C. sakazakii strains and the selected LAB strains, which inhibited C. sakazakii and did not exhibit atypical antibiotic resistance, were grown in Luria-Bertani or de Man Rogosa Sharpe broth, respectively, containing 1% dextrose or 1% commercial prebiotic (w/v) to compare their ability to metabolize the prebiotic product. Overall, based on the growth inhibition of C. sakazakii, antibiotic susceptibility, and prebiotic metabolism, 6 of the 16 LAB were chosen to be part of a potential synbiotic. This study has provided valuable information that will help with the development of a synbiotic that can be used in powdered infant formula to reduce the potential for C. sakazakii-related illnesses in infants.

Application of a High-Throughput Targeted Sequence AmpliSeq Procedure to Assess the Presence and Variants of Virulence Genes in Salmonella.

We have developed a targeted, amplicon-based next-generation sequencing method to detect and analyze 227 virulence genes (VG) of Salmonella (AmpliSeqSalm_227VG) for assessing the pathogenicity potential of Salmonella. The procedure was developed using 80 reference genomes representing 75 epidemiologically-relevant serovars associated with human salmonellosis. We applied the AmpliSeqSalm_227VG assay to (a) 35 previously characterized field strains of Salmonella consisting of serovars commonly incriminated in foodborne illnesses and (b) 34 Salmonella strains with undisclosed serological or virulence attributes, and were able to divide Salmonella VGs into two groups: core VGs and variable VGs. The commonest serovars causing foodborne illnesses such as Enteritidis, Typhimurium, Heidelberg and Newport had a high number of VGs (217-227). In contrast, serovars of subspecies not commonly associated with human illnesses, such as houtenae, arizonae and salame, tended to have fewer VGs (177-195). Variable VGs were not only infrequent but, when present, displayed considerable sequence variation: safC, sseL, sseD, sseE, ssaK and stdB showed the highest variation and were linked to strain pathogenicity. In a chicken infection model, VGs belonging to rfb and sse operons showed differences and were linked with pathogenicity. The high-throughput, targeted NGS-based AmpliSeqSalm_227VG procedure provided previously unknown information about variation in select virulence genes that can now be applied to a much larger population of Salmonella for evaluating pathogenicity of various serovars of Salmonella and for risk assessment of foodborne salmonellosis.

Salmonella enterica subsp. enterica virulence potential can be linked to higher survival within a dynamic in vitro human gastrointestinal model.

Salmonella enterica subsp. enterica is one of the leading causes of human foodborne infections and several outbreaks are now associated with the consumption of fresh fruit and vegetables. This study aims at evaluating whether Salmonella virulence can be linked to an enhanced ability to survive successive digestive environments. Thirteen S. enterica strains were selected according to high and low virulence phenotypes. Lettuce inoculated separately with each S. enterica strain was used as food matrix in the TNO gastrointestinal model (TIM-1) of the human upper gastrointestinal tract. During the passage in the stomach, counts determined using PMA-qPCR were 2-5 logs higher than the cultivable counts for all strains indicating the presence of viable but non-cultivable cells. Bacterial growth was observed in the duodenum compartment after 180 min for all but one strain and growth continued into the ileal compartment. After passage through the simulated gastrointestinal tract, both virulent and avirulent S. enterica strains survived but high virulence strains had a significantly (p = 0.004) better average survival rate (1003 %-3753 %) than low virulence strains (from 25 % to 3730%). The survival rates of S. enterica strains could be linked to the presence of genes associated with acid and bile resistance and their predicted products. The presence of single nucleotide polymorphisms may also impact the function of virulence associated genes and play a role in the resulting phenotype. These data provide an understanding of the relationship between measured virulence potential and survival of S. enterica during dynamic simulated gastrointestinal transit.

Genomic and phenotypic analysis of SspH1 identifies a new Salmonella effector, SspH3.

Salmonella is a major foodborne pathogen and is responsible for a range of diseases. Not all Salmonella contributes to severe health outcomes as there is a large degree of genetic heterogeneity among the 2,600 serovars within the genus. This variability across Salmonella serovars is linked to numerous genetic elements that dictate virulence. While several genetic elements encode virulence factors with well-documented contributions to pathogenesis, many genetic elements implicated in Salmonella virulence remain uncharacterized. Many pathogens encode a family of E3 ubiquitin ligases that are delivered into the cells that they infect using a Type 3 Secretion System (T3SS). These effectors, known as NEL-domain E3s, were first characterized in Salmonella. Most Salmonella encodes the NEL-effectors sspH2 and slrP, whereas only a subset of Salmonella encodes sspH1. SspH1 has been shown to ubiquitinate the mammalian protein kinase PKN1, which has been reported to negatively regulate the pro-survival program Akt. We discovered that SspH1 mediates the degradation of PKN1 during infection of a macrophage cell line but that this degradation does not impact Akt signaling. Genomic analysis of a large collection of Salmonella genomes identified a putative new gene, sspH3, with homology to sspH1. SspH3 is a novel NEL-domain effector.